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ABSTRACT 

We investigate how the hierarchical merging of dark matter halos, the radiative cooling of baryons, and the 
energy feedback from supernovae and active galactic nuclei or quasars combine to govern the amount and the 
thermal state of the hot plasma pervading groups and clusters of galaxies. We show that, by itself, supernova 
preheating of the external gas flowing into clusters falls short of explaining the observed X-ray scaling relations 
of the plasma luminosity Lx or the plasma entropy K versus the X-ray temperature T. To account for the scaling 
laws from rich to poor clusters takes preheating enhanced by the energy input from active galactic nuclei. In 
groups, on the other hand, the internal impacts of powerful quasars going off in member galaxies can blow 
some plasma out of the structure. So they depress Lx and raise K to the observed average levels; meanwhile, 
the sporadic nature of such impulsive events generates the intrinsic component of the wide scatter apparent in 
the data. The same quasar feedback gives rise in groups to entropy profiles as steep as observed, a feature hard 
to explain with simple preheating schemes. Finally, we argue a close connection of the Lx - T or the K—T 
relation with the M. - a correlation between the host velocity dispersion and the masses of the black holes, 
relics of the quasar activity. 

Subject headings: galaxies: clusters: general - quasars: general - X-rays: galaxies: clusters 



1. INTRODUCTION 

The hot medium pervading many single galaxies and 
most groups and clusters shines in X rays by thermal 
bremsstrahlung and line emission, see Sarazin (1988). Simple 
conditions are found to prevail in rich clusters. 

These emit huge powers Lx oc n 1 \fT R 3 ~ 10 44 - 10 45 erg 
s" 1 in X rays; the temperatures kT w 5 keV, measured from 
the continuum and from high excitation lines, are close to 
the virial values kT v w GMm p /lQR in the gravitational wells 
mainly provided by dark matter (DM) masses M w 10 15 M Q 
within sizes R of a few Mpcs 1 . The inferred gas number den- 
sities decline outward from central values n w 10~ 3 cm" 3 ; so 
this medium with low n and high T satisfying kT /e 2 n x ^ ~ 
10 12 constitutes the best ion-electron plasma in the Universe 
ever, the intracluster plasma or ICP 

Such a medium is apparently simple on the following ac- 
counts. Microscopically, it is constituted by pointlike par- 
ticles in thermal equilibrium. At the macroscopic end, the 
overall baryonic fractions resulting from the ICP densities and 
radial distributions inventoried in many clusters (White et al. 
1993, Allen & Fabian 1994) come to values m/M w 0. 16; this 
is close to the cosmic ratio f2j / ft M of baryons to DM obtained 
in the current Concordance Cosmology from the parameters 
Q, b w 0.044 and fl M » 0.23 (see Bennett et al. 2003) 2 . In ad- 
dition, the chemical composition is reasonably constant from 
cluster to cluster, and close to one-third of the solar value (see 
Matteucci 2003). 

However, surprises arise in moving from rich clusters to- 
ward poor groups. In fact, similarly simple conditions holding 

1 Throughout the paper k is the Boltzmann constant, G is the gravitational 
constant, m p is the proton mass, and e is the electron charge. 

2 This will be used throughout the paper, with its additional parameters 
£2a ~ 0.73 for the dark energy density, and Hq ~ 70h7o km s _I Mpc - ' for 
the Hubble constant. 



in the intragroup plasma (IGP) would imply that the luminosi- 
ties retain the gravitational scaling Lx oc T 2 (Kaiser 1986). 
This would apply if the IGP passively followed the DM evo- 
lution, and retained the key cluster behaviors: m/M ~ const, 
i.e., densities n proportional to the DM mass density p; and 
temperatures T close to the virial value T v ocM/Roc M 2 / 3 p 1 ' 3 . 

Instead, the luminosities recently detected or revised 
(Horner 2001; O' Sullivan, Ponman & Collins 2003; Osmond 
& Ponman 2004) are lower by factors of 10 _1 to 10~ 2 , see 
Fig. 1 . The figure also shows how the emissions from poor 
groups and large galaxies scatter widely and often downward, 
a feature of a largely intrinsic nature (Mushotzky 2004). 

So in such smaller structures the plasma is surprisingly un- 
derluminous and hence underdense. This is an even more 
surprising result, considering that in the standard hierarchi- 
cal cosmogony (see Peebles 1993) such earlier condensations 
ought to be denser, if anything. Moreover, for kT < 2 keV 
the pinch of highly excited metals contributes important line 
emissions that imply a flatter Lx oc T, if anything. How the 
observed steep decrease may come about constitutes a widely 
debated issue. 

Our proposal centers on the energy gained or lost by the 
baryons through several processes: the gravitational heating 
driven by the merging events that punctuate the hierarchical 
growth of DM condensations ("halos", see Peebles 1993); the 
radiative cooling of the baryons; and the energy fed back to 
baryons when they partly condense within galaxies into mas- 
sive stars then exploding as Type II supernovae (SNe), or ac- 
crete onto a central supermassive black hole (BH) energizing 
an active galactic nucleus (AGN) or a quasar. 

The paper is organized as follows. In § 2 we use a telling 
quantity, the plasma entropy, to show that energy feedback 
from astrophysical sources is needed to explain its high levels 
in poor clusters and groups. In § 3 we show that preheating 
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by SNe alone is not enough. In § 4 we consider preheating 
enhanced by AGNs, compute the resulting X-ray scaling re- 
lations, and critically discuss the outcomes of this approach. 
In § 5 we also consider the internal impacts from quasars; 
we compute their effect on the X-ray scaling relations and on 
the entropy profiles in groups, and show how these solve the 
shortcomings of all external preheating scenarios. In § 6 we 
highlight and discuss the main features of our comprehensive 
approach. 

Auxiliary computations and derivations are given in the Ap- 
pendices. In Appendix A we derive Eq. (2) of the main text 
from the hydrostatic equilibrium, and give handy approxima- 
tions. In Appendix B we reformulate the classic Rankine- 
Hugoniot jump conditions in a general form that is also valid 
for accretion shocks/layers, and derive Eqs. (5) and (6) of the 
main text. In Appendix C we develop a new family of self- 
similar hydrodynamic solutions describing the blast waves 
driven by the internal impacts of quasars, and extensively used 
in §5. ' 

2. X-RAY LUMINOSITY AND ENTROPY 

This paper deals with the energy budget of the baryons. The 
latter experience gravitational heating to T s» T v (discussed in 
detail in § 4) as they fall into the hierarchically growing DM 
structures. Nongravitational energy losses or inputs deplete 
the baryon density; this is because losses trigger baryon con- 
densation to stars, while inputs cause outflow from and hinder 
inflow into newly forming structures. 

All such processes are probed with the adiabat K = 
kTn~ 2 ' 3 oc e 2s / 3/; , a direct measure of the specific entropy s 
(see Bower 1997; Balogh, Babul & Patton 1999). The levels 
of the adiabat K are linked to Lx by the inverse relation 



Ko,L x l/3 T^ 



(1) 



which obtains at z ~ on eliminating n between their respec- 
tive expressions; note that T 5 / 3 goes over to T 4 ^ for impor- 
tant line emission, and that we are neglecting the weakly T- 

1/3 

dependent shape factors for K and L x . 

Clearly, K stays constant under adiabatic transformations of 
the plasma. Gravitational heating would set the scaling K oc T 
(corresponding to Lx oc T 2 after Eq. [1]), but Fig. 2 shows the 
data for decreasing T to deviate substantially upwards; this 
indicates that additional nongravitational processes occurred 
during a structure's merging history. The present paper in- 
vestigates how these processes affect the adiabat K = K 2 n(r), 
namely, the level K 2 at the virial radius r = R, and the inner 
profile K(r). 

The relation of these quantities to the density run n(r) of the 
plasma in hydrostatic equilibrium is derived in Appendix A, 
and reads 



n(f) 

«2 



— = «" 3 / 5 (r) 



3/2 



(2) 



the boundary condition at r = R is given by n 2 = (kT 2 /K 2 fl 2 , 
but K 2 will be related to T 2 in § 4. The parameter j3 = T v /T 2 is 
the ratio of the DM to the thermal plasma scale height in the 
gravitational potential 4>(r); the latter is normalized to the one- 
dimensional dispersion a 2 = kT v //j,m p , with fi w 0.6 for the 
nearly cosmic composition of the plasma (Cavaliere & Fusco- 
Femiano 1976). For a polytropic entropy distribution k(t) oc 
n(r) r-5 ' 3 with uniform index T = 5/3 + d\nn/d\nn, Eq. (2) 
yields the familiar form n = n 2 [1 + (T - I) (3 A(j>/T] l ^ r - l) in 
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FIG. 1. — Integrated X-ray luminosity Lx vs. X-ray temperature T. Data 
for clusters (crosses) are from Horner (2001), for groups (circles) from Os- 
mond & Ponman (2004), and for early-type galaxies (stars) from O' Sullivan 
et al. (2003). The dotted line represents the gravitational scaling, with line- 
emission included. The strip (with 2 a width provided by the merging histo- 
ries) illustrates our results for SN preheating with IcAT = 1 keV per particle, 
as discussed in § 3. 



potential drop Aqb inward of R; the isothermal 
•p A<£ f or y = 1 provides the standard model to fit 



terms of the 

limit n = «2r for V = 1 provides 1 
the X-ray surface brightness profiles, which yields values (3 ~ 
0.7 in rich clusters. In particular, we use Eq. (2) to compute 
integrated luminosities L x oc j drr 2 n 2 {r)T l l 2 {r) and central 
entropies Kq,\ at r = 0.1 R, that we compare with the data. 

Throughout the paper we will make use of "semianalytic" 
techniques. This is because the nongravitational processes af- 
fecting K include energy inputs and radiative losses, which 
interplay in complex patterns with gravitational heating; so 
even numerical simulations based on advanced A^-body and 
hydrocodes are driven to, or beyond their present limits, 
and have to borrow from semianalytic models much subgrid 
physics (see discussion by Borgani et al. 2002). 

Concerning radiative losses, they do operate within galax- 
ies to remove low-entropy gas by condensing it into stars, a 
process that we include in our semianalytic modeling. But 
extensive cooling as needed to substantially raise the resid- 
ual ICP/IGP entropy or depress Lx would produce too many, 
unseen stars (see Voit & Bryan 2001; Sanderson & Ponman 
2003). On the other hand, cooling triggers catastrophic in- 
stabilities unless closely restrained by other processes feed- 
ing energy back to baryons (see Blanchard, Valls-Gabaud & 
Mamon 1992); so energy additions AE > are mandatory 
anyway, and will constitute our focus next. 

3. THE NEED FOR AGN FEEDBACK 

Obvious energy sources are provided by Type II SN ex- 
plosions; do they contribute enough energy feedback? SNe 
provide energies E$n ~ 10 51 ergs with an occurrence t]sn ;$ 
5 x 10~ 3 per solar mass condensed into stars, the latter value 
being calibrated so as to include the yield of strong winds 
from young hot stars (Bressan, Chiosi & Fagotto 1994). 

Such outputs may be coupled to the surrounding gas at lev- 
els fsN ;$ 5 when cooperative SN remnants propagation takes 
place, as in the case of starbursts, to drive subsonic galactic 
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FIG. 2. — Central entropy ATo. i (at r 7a 0. 1 if ) vs. X-ray temperature T. Data 
for clusters and groups are from Ponman et al. (2003): circles mark individual 
systems and squares refer to binned data. The dotted line represents pure 
gravitational heating. The strip (with 2 a width provided by the merging 
histories) illustrates our results for SN preheating with kAT = | keV per 
particle, as discussed in § 3. 



winds (see Matteucci 2003). Then the integrated thermal in- 
put attains the maximal level (Cavaliere, Lapi & Menci 2002) 

kAT = f SN E SN Vsn—<^ keV particle" 1 (3) 

3 m 4 

in groups with stellar to gas mass ratios up to m*/ m s=s i (see 
David 1997); somewhat smaller values obtain in clusters (e.g., 
Lin, Mohr & Stanford 2003). 

The SN feedback and the originating star formation are de- 
scribed by semianalytic models, in particular that of Menci & 
Cavaliere (2000). They base their model on the DM merging 
histories, i.e., the hierarchical buildup of a galaxy or a group 
with their DM and baryonic contents through merging events 
with comparable or smaller partners, down to nearly smooth 
inflow (Lacey & Cole 1993). The model, in addition, specifies 
how the baryons are cycled between the cool, stellar, and hot 
phases; the latter contributes to the ICP/IGP, while the former 
two phases yield the stellar observables. 

It is found that most of a structure's DM mass M (and of 
the IGP mass likewise) is contributed to its main progenitor 
by merging partners with masses M' < M/3 and related virial 
temperatures TJf < 0.6 T v (Cavaliere, Menci & Tozzi 1999). 
The smaller lumps have shallower gravitational wells and pro- 
duce more star-related energy on scales closer to their dynam- 
ical times; so they are more effective in heating up their gas 
share to temperatures Tj + AT. During each subsequent step 
of the hierarchy forming larger groups or clusters, the exter- 
nally preheated gas (see Muanwong et al. 2002) will be hin- 
dered from flowing in and contributing to the IGP or ICP So 
under any model depleted densities will be propagated some 
steps up the hierarchy. 

In sum, SNe make optimal use of their energy in preheat- 
ing the IGP. However, their input kAT w | keV per particle 
turns out to cause only limited luminosity depressions or en- 
tropy enhancements, as shown by the light strips in Figs. 1 
and 2. The result may be understood by referring to the simple 
isothermal case where L x oc n\J drr 2 e 2i:sA ^ applies; in mov- 



ing from rich to poor clusters n decreases, governed mainly 
by the decreasing exponential e 2f3A ^ (as visualized by Fig. 5 
of Cavaliere & Lapi 2005). But the normalized DM potential 
A(f> deepens because of the increased concentration (an intrin- 
sic feature discussed in Appendix A). To offset this trend and 
provide constant or decreasing density, it is clearly required 
that f3 w T V /(T V + AT) be lowered from the cluster value by 
a sufficiently strong preheating AT; in detail, the approxima- 
tion (3 i=a 0.7- AT /T holds; see Appendix B and in particu- 
lar Eq. (B9) with 9 = T/AT. Numerically, the requirement 
comes to kAT > 0.5 keV per particle for any significant lu- 
minosity depression in a poor cluster; yet more preheating is 
required with polytropic plasma distributions. On the other 
hand, including Type la SNe still does not meet the above 
requirement (see Pipino et al. 2003). 

In view of these SN limitations, in the rest of the paper we 
concentrate on the stronger feedback provided by quasars and 
active galactic nuclei (see Valageas & Silk 1999; Wu, Fabian 
& Nulsen 2000; Yamada & Fujita 2001; Nath & Roychowd- 
hury 2002). These sources are kindled when sizeable amounts 
of galactic gas, triggered by mergers or interactions of the host 
with companion galaxies, are funneled downward from kpc 
scales; they not only form circumnuclear starbursts but even- 
tually trickle farther down to the very nucleus (see Menci et 
al. 2004) and accrete onto a central supermassive BH. 

4. EXTERNAL PREHEATING FROM AGNs 

On accreting the BH mass M,, the integrated energy input 
to the surrounding plasma comes to values 

kAT^fnc 2 ^* I keV particle-, (4) 
3 4Mb m 2 

easily larger than for SNe. We have used the standard mass- 
energy conversion efficiency (/« 10 _1 , and the locally ob- 
served ratio M,/Mb w 2 x 10" 3 ofBH to galactic bulge masses 
(Merritt & Ferrarese 2001); the factor \ accounts for the bulge 
mass observed in blue light compared to that integrated over 
the star formation history (Fabian 2004a). Finally, we will 
adopt values /w5x 10~ 2 for the fractional AGN output actu- 
ally coupled to the surrounding gas, on the grounds discussed 
next. 

The 10% radio-loud AGNs directly produce considerable 
kinetic or thermal energies in the form of bubbles and jets (see 
Forman et al. 2004), but statistics and nonspherical geometry 
reduce their average contribution to /. On the other hand, in 
the 90% radio-quiet AGNs a small coupling is expected for 
the radiative output because of the flat spectrum and low pho- 
ton momenta; the observations of wind speeds up to v w w 0.4c 
suggest values around v w /2c w 10 _1 associated with covering 
factors of order 10" 1 (see Chartas, Brandt & Gallagher 2003; 
Pounds et al. 2003). We shall see that average values around 
/ w 5 x 10" 2 are consistent not only with the X-ray observa- 
tions of the IGP, but also with the mainly optical observations 
of the relic BHs in galaxies. 

Considering that the AGN activity closely parallels the star 
formation in spheroids (Franceschini et al. 1999; Granato et 
al. 2004; Umemura 2004), we add the AGN energy injections 
to SN's to obtain preheating energies up to kAT rts | keV 
per particle. Such a combined value produces a sizeable step 
toward the locus of the data, as shown in Figs. 3 and 4 by the 
heavy strips. How we obtain these is explained next. 

Toward this purpose it is convenient to discuss further the 
modus operandi of the external preheating. During the for- 
mation of a DM structure, outer lumps and the associated 
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gas flow in together; but just inside R the smaller and/or less 
bound gas bunches are promptly stripped away from their DM 
hosts, while gaining entropy (Tormen, Moscardini & Yoshida 
2004, see their Fig. 8). The outcome constitutes a complex 
patchwork of shocks of all sizes, comprised within an outer 
layer with thickness 6 < 10 _1 R wherein most of the entropy 
rise takes place. 

The net result is close to that computed from considering 
a coherent accretion shock, roughly spherical and located at 
r ~ R as considered by Cavaliere et al. (1999), Dos Santos 
& Dore (2002), and Voit et al. (2003). In fact, across the 
layer we may retrace the classic Rankine-Hugoniot derivation 
based on the conservation of mass, momentum and energy 
for plasma particles with 3 degrees of freedom; for a reason- 
ably thin layer we recover the standard entropy jump across a 
shock 



K 2 =K l 5/3 



2(6-1)+ d 4 (6 



■if + t 



-2/3 



(5) 



to within 0(5/ R) accuracy, independently of geometrical de- 
tails (as discussed in Appendix B). To within the same accu- 
racy, the strength parameter 9 = T 2 /T l7 i.e., the ratio of the 
down- to the upstream temperature is linked by 



o= 5 -^ + 7 



16 



8 16M 2 



(6) 



to the Mach number M = (3 pm p v\/5kT\) l l 2 of the flow ve- 
locity V\ relative to the shock/layer. Henceforth we refer to 
"shocks", but we also include thin layers. 

We now consider in closer detail the combined preheating 
by SNe and AGNs; this comes into play through Ti=Tj + AT 
that enters K\ and 6 in Eqs. (5) and (6). We average the result- 
ing K 2 over the full structure's merging history that includes 
the distributions of progenitor masses M' or related T',; to this 
purpose we implement, as in Cavaliere et al. (1999), the con- 
ditional probabilities and the merging rates from the standard 
cold DM cosmogony as given by Lacey & Cole (1993). This 
straightforward if laborious procedure (which is dominated 
by the smaller partners and so further validates Eqs. [5] and 
[6], as discussed in Appendix B) is made semianalytically, 
and yields the heavy strips in Figs. 3 and 4. Their width il- 
lustrates the variance (at 95% probability level) around the 
mean value, induced mainly by the merging stochasticity; the 
smooth, low-power AGN activity considered here does not 
contribute much additional scatter. 

Next we discuss why the results fit clusters better than 
groups. In very rich clusters the infall velocities v\ w 
2.1 (kT v / /imp) 1 / 2 are large, and v\ ~ 4vi/3 is larger yet (see 
Appendix B). These velocities are dominantly supersonic, ex- 
cept for the few major lumps which carry warm gas deep into 
the structure (as observed by Mazzotta et al. 2002) and con- 
tribute little to prompt entropy gains. So the effective shocks 
are uniformly strong with ~ /j,m p v 2 /3 kT\ 3> 1 , see Eq. (B7). 
Such conditions in Eq. (5) yield K2/K1 ~ 9/4 2 / 3 , correspond- 
ing to nearly constant «2/«i ~ 4, see Eq. (B7); they also yield 
a nearly constant value of (3 = T v /0 T\ ~ 3 kT v / pjn p v\ ~ | . In 
other words, here we find pure gravitational heating at work 
to enforce K oc T or L x oc T 2 . 

The related, raising entropy profiles reflect the history of 
progressive depositions of shells undergoing stronger and 
stronger shocks during the hierarchical growth to a rich clus- 
ter; in fact, in the outer regions we find (tar 1,1 and T w 1.1. 
We derive these values by reducing to bare bones the model 
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FIG. 3. — Integrated X-ray luminosity Lx vs. X-ray temperature T. Data, 
dotted line and light shaded strip as in Fig. 1 . The heavy shaded strip (with 2 
a width provided by the merging histories) illustrates our results for external 
preheating when including the AGN contribution to a total kAT = | keV per 
particle, as discussed in § 4. Our results for the internal impacts from quasars 
are illustrated by the solid (ejection model) and dashed (outflow model) lines, 
see § 5 for details. The coupling level of the quasar output to the ambient 
medium is / = 5 X 10 -2 . 



of Tozzi & Norman (2001). We adopt in full the Concordance 
Cosmology, and nearly self-similar hierarchical collapse with 
a constant perturbation power index around —1.2 appropriate 
for rich clusters (see Padmanabhan 2003); these conditions 
imply the last accreted shell to add a mass AM oc M on top 
of the mass M oc (1 +z)~ 3 2 virialized at z > \ (see also Lapi 
2004). In the process, the entropy K oc T /n 2 ! 3 grows because 
the strong shocks prevailing in rich clusters yield not only 
T w T v and «2 — 4n\ as above, but also n\ oc p oc (1 +z) 2 9 con- 
sidering the appropriate collapse threshold; in terms of m oc M 
this translates into k oc m 2 ! 3 (1 + z) _1 oc m. 

For the ICP in equilibrium, the radial entropy profile n(r) 
corresponding to this distribution n(m) is found as follows. 
In the outskirts we approximate the entropy profile as n(r) oc 
r a , with a to be determined; then Eq. (A3) implies n(r) oc 
r -3( Q +2)/5^ an( j we b tam K ( r ) oc m(< r ) 5a / 3 ( 3 -") n consid- 
ering that m(< r) = 4nm p J dx^n(x) holds. Requiring con- 
sistency with the entropy distribution n oc m derived above, 
we obtain a w 1.1. So our final results read «(r) oc r 11 and 
n(r) oc r" L9 , which accord with the data by Ponman, Sander- 
son & Finoguenov (2003), and with the simulations by Tor- 
natore et al. (2003); the related value V w 1 . 1 agrees with that 
observed by Ettori & Fabian (1999) and De Grandi & Molendi 
(2002). 

In poor clusters, on the other hand, the infall is slower 
with v\ ~ 4vi/3 + 5kTi/4pm p v\, see Appendix B. Now the 
inflow is less supersonic, and the accretion shocks are eas- 
ily modulated by the preheating temperature to a strength 
6 ~ nm p v 2 /3kT\ +3/2, see Eq. (B8). Less entropy is pro- 
duced by these intermediate shocks, while an additional con- 
tribution is just carried in with the warm inflowing gas, to 
yield K2/K1 ~ (6* + 5/8)/4 2 / 3 ; correspondingly, the boundary 
densities are lowered to «2/"i — 4- 15/40, see Eq. (B8). In 
addition, the density profiles are now just flatter than in rich 
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FIG. 4. — Central entropy Kqa vs. X-ray temperature T. Data, dotted 
line and light shaded strip as in Fig. 2. The heavy shaded strip (with 2 a 
width provided by the merging histories) illustrates our results for external 
preheating when including the AGN contribution to a total kAT = | keV per 
particle, as discussed in § 4. Our results for the internal impacts from quasars 
are illustrated by the solid (ejection model) and dashed (outflow model) lines, 
see § 5 for details. As before, the coupling level is / = 5 X 1(T 2 . 



clusters, since n(r) jm is appreciably decreased with preheat- 
ing levels kAT w | keV per particle that satisfy the condition 
derived at the end of § 3; these are effective in lowering all 
densities, hence in depressing L x and enhancing K. 

In groups and galaxies these preheating levels are enough to 
cause smoother, transonic inflows and weak shocks with 8 ~ 
1, yielding small jumps K2/K1 ~ 1 + 5 (6 1 — l) 3 / 6 and «2/«i — 
1, see Eq. (BIO). The X-ray scaling relations produced by 
the combined external preheating of SNe and AGNs (heavy 
strips in Figs. 4 and 3) are in marginal agreement with the 
trends in the data, while the wide scatter is still unaccounted 
for. Moreover, the weak shocks so produced would imply 
nearly flat profiles K,(r), which often are not observed (Pratt & 
Arnaud 2003; Rasmussen & Ponman 2004). 

The problem with such isentropic profiles would be aggra- 
vated and propagated to poor clusters while the entropy would 
be raised too much, if one tuned high the AGN preheating, 
much above the level \ keV per particle given by Eq. (4). 
From the previous relations the problem is easily seen to de- 
velop even before solving the luminosity issue in groups. 

5. INTERNAL IMPACTS OF QUASARS 

But right in groups and galaxies the impulsive inputs by 
powerful quasars take over, providing from inside an addi- 
tional impact on the IGP that can cause outflow or ejection. 
For this to occur, two energies compete: the overall input 
A£k 2 x 10 62 / (M./10 9 M©)(1 +zT 3/2 ergs provided by a 
quasar on accreting the mass M, within the host dynamical 
time td w 10 8 yr set by mergers or interactions; the (absolute) 
total energy E w2x 10 61 (fcr/keV) 5 / 2 (1 +z)" 3//2 ergs residing 
in the equilibrium IGP (Lapi, Cavaliere, & De Zotti 2003). 

The relevant ratio 



is small in clusters but increases toward groups, and ap- 
proaches unity in poor groups with kT w 1 keV to attain a 
few in large galaxies with kT w \ keV. Within the central 
kpc of such structures the quasar launches a piston (see King 
2003; Granato et al. 2004) that drives through the surround- 
ing plasma a blast wave bounded by a leading shock at r = R s 
(see Fig. 5). These blasts constitute effective, quasi-isotropic 
means to propagate energy far away from the central source. 

While the latter shines, the blast affects the plasma out to 
the distance R s where the initial energy E(< R s ) is compa- 
rable to the cumulative input AE(t). In fact, the condition 
AE(t) /E(< R s ) = const defines the self-similar propagation of 
the blast and the motion R s (t) of the leading shock; as shown 
in Appendix C, the result is R s (t) oc t 2// " in an initial distribu- 
tion n{r) oc r _a) (2 < ui < 2.5) for the plasma density under the 
energy input AE(t) oc f 2 ( 5 - 2w )/ w The simplest flow obtains 
with u = 2 (the standard isothermal sphere) implying a source 
power L(t ) = d AE jdt = const; when u> > 2 applies the power 
L(t ) oc t 5 ( 2-w V w declines, a useful means to describe the quasar 
fading because of its own feedback on the accreting gas. 

This new family of hydrodynamic solutions, proposed by 
Cavaliere et al. (2002), is described in Appendix C; it is used 
below to evaluate mass loss from, and entropy distribution 
into, the structures. The solutions include the restraints set to 
gasdynamics by a nonzero initial pressure p(f) oc r 2(1_a,) and 
by the DM gravity; so not only do they imply a well-defined 
E(< R s ), but they also cover the full range of blast strengths 
from weak in clusters to strong in galaxies, depending on the 
magnitude of the key parameter AE/E. 

As AE/E increases, so does the Mach number M. of the 
leading shock; their relation is shown in Fig. C2 and is ap- 
proximated by M 2 ~ 1 + (AE /E) for AE/E < 2 in the simple 
u> = 2 model. Meanwhile, the ratio of the kinetic to the thermal 
energy ranges up to 2, see Eq. (C17). Correspondingly, within 
the time tj increasing plasma amounts are driven beyond the 
virial radius R of a large galaxy or a poor group; in the simple 
u> = 2 model the fractional mass ejected or flowed out is well 
approximated by Am/m ~ AE/2E, see also Table 1. 



n(r) 



\ L 



n. 



A£ / M. /kT_ 

E ~ ' 5 x 10" 2 1O 9 M VkeV 



-5/2 



(7) 



R p R s R 

FIG. 5. — An outline of the density distribution during the propagation of a 
quasar-driven blast throughout the equilibrium plasma. By the dashed line we 
represent the initial density run n oc r in the volume already evacuated by 
the blast, by the thick solid line the perturbed density in the blast (specified in 
Fig. CI), and by the thin solid line the still unperturbed density. The perturbed 
flow is confined between the trailing piston at R p and the leading shock at R s ; 
in our illustration this is still far from the virial radius R. 
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TABLE 1 

Relevant Blast Wave Quantities (u> = 2) 



AE/E 


M 


(P)/Pi 


1-Am/m 




0.3 


1.2 


3.6 


0.92 


0.94 


1 


1.5 


4.6 


0.58 


0.86 


3 


1.9 


6.3 


~0 





These results turn out to be nearly independent of the spe- 
cific mode for mass loss; in particular, we compare two ex- 
treme cases. In the first one ("ejection"), we take Am to be 
the mass in the blast driven outside of the virial radius R at 
t = td by the blast kinetic energy. In the second case ("out- 
flow"), we adopt constant pressure as boundary condition at 
r = R, to obtain new densities n' oc (T+ AT)" 1 ; now Am is 
the mass flowed out of the structure due to the extrathermal 
energy Ar deposited by the blast. Beyond model details, we 
find that the mass loss closely obeys Am/m ~ AE/2E. 

In both cases, after the passage of the blast the IGP will re- 
cover hydrostatic equilibrium, described by Eq. (2). But all 
new densities n' will be depleted by the factor 1 - Am/m be- 
low the initial value already affected by the preheating from 
SNe and AGNs. In addition, in both cases the extrather- 
mal energy deposited by the blast lowers the values of 0'; 
this is given in Table 1, and may be understood in terms of 
0' 10 w T/(T + AT). The resulting L x oc (1-Am/m) 2 , in- 
cluding the appropriate T-dependent shape factor, is shown 
by the solid and dashed lines in Fig. 3 for both models, ejec- 
tion and outflow, discussed above. 

The IGP entropy is increased by quasar-driven blasts. 
While these sweep through the plasma, a moderate produc- 
tion takes place across the leading shock, and leads to a jump 
K2/K1 given again by the general Eqs. (5) and (6). In the equi- 
librium recovered after the plasma mass loss Am/m caused by 
ejection or outflow, the entropy is further enhanced to read 

K!,/K2= (l-Am/m)" 2/3 . (8) 

For example, in a group with kJ = | keV the combined pre- 
heating by SNe and AGNs yields entropy levels correspond- 
ing to 100 keV cm 2 . This is raised to 180 keV cm 2 by 
the internal blast driven by a quasar deriving from a BH of 
10 9 M© with coupling / w 5 x 10" 2 , that produce AE/E w 1. 
The resulting central entropy Kq,\ including the appropriate 
T-dependent shape factor is shown by the solid and dashed 
lines in Fig. 4 for both models, ejection and outflow, discussed 
above. 

Relatedly, the entropy profiles are steep after the blast pas- 
sage. They arise when the entropy produced in the blast, 
clearly piled up toward the leading shock, is redistributed in 
the recovered equilibrium. We find 

K(m) oc m 4//3 (9) 

to hold in terms of the plasma mass m swept up by the blast 
(see Fig. CI and Eq. [Cll]). In the adiabatically recovered 
equilibrium we require the entropy distribution n'(m) to equal 
n(m) oc m 4 / 3 , and proceed in analogy with the technical steps 
used in § 3. 

In detail, let us approximate k' oc r a , with a to be de- 
termined; then Eq. (A3) implies n' oc r - 3 < a+2 )/ 5 ? and thus 
k' oc m W 3 ( 3 -") follows. Requiring this to be consistent with 



the entropy distribution k oc m 4//3 in the blast (as anticipated 
above) yields a s=s 1.3. In other words, the blast acting from 
inside leaves in the readjusted plasma a strong imprint of 
its own entropy distribution, in the form of a steep profile 
n\r) oc r 1 ' 3 consistent with the data. 

6. DISCUSSION AND CONCLUSIONS 

In this paper we have used pilot semianalytic modeling 
to show that the energy fed back to baryons by AGNs and 
quasars is essential to fit the recent X-ray data. We find that 
the AGN external preheating dominates over SN's to yield the 
scaling laws Lx oc T 3 and K oc T 2 / 3 related by Eq. (1), that 
constitute fitting trends in clusters if still marginal for groups. 
But in groups and in large galaxies the quasar impulsive feed- 
back acting from inside takes over, to eject some plasma and 
further depress Lx, while enhancing K and originating non- 
isentropic profiles. 

We stress that the energy added by AGNs plays an inverse 
role in preheating and in ejection/outflow. This is because 
in moving from clusters to groups the ratio AE/E of the 
added vs. the equilibrium energy is bound to increase. This 
causes relatively higher external preheating, warmer inflows 
and weaker accretion shocks; on the other hand, it drives 
stronger internal blasts causing more mass ejection/outflow. 
In parallel, the leading shocks of the blasts replace the accre- 
tion shocks in the role of increasing the outer entropy produc- 
tion, so originating comparably steep entropy profiles. 

Such internal effects, however, are to saturate in galaxies 
because the values of AE/E there are limited, lest the im- 
pulsive quasar feedback ejects so much of the surrounding 
baryons as to stop the BH accretion altogether (see Silk & 
Rees 1998); the saturation is what in our calculation yields 
the lower elbows of the solid and dashed lines presented in 
Fig. 3. As argued in Cavaliere et al. (2002), the system consti- 
tuted by the BH and the surrounding baryons self-regulates to 
the verge of unbinding; the condition AE w E directly yields 
M,«5x 10 8 M q (//5 x lO^r^cr/SOOkms" 1 ) 5 in terms of 
the DM velocity dispersion a. In turn, the latter is found to 
correlate less than linearly with the velocity dispersion cr* of 
the host galactic bulge (Ferrarese 2002; Baes et al. 2003; 
A. Pizzella et al. 2005, in preparation); so the BH mass ap- 
proaches M. OC (J 4 . 

Thus for the same value / w 5 x 10~ 2 indicated by the av- 
erage X-ray data for groups we agree (within the observed 
scatter, see Tremaine et al. 2002) with the observations of the 
uppermost relic BH masses in the bulges of many local and 
currently inactive galaxies. Similar BH masses may be also 
contributed by an initial, supercritical accretion phase (as dis- 
cussed by King 2003) launching the piston that in turn drives 
the far-reaching blasts described above. Subsequently, our 
outgoing self-similar blasts with AE(t)/E(< R s ) w const stay 
tuned to the condition AE w E. Specifically, for uj — > 2.5 
not only E{< R s ) oc R 5 ~ 2uJ — > const holds but also AE oc 
? 2(5-2w)/w _^ const applies, consistent with fading quasar out- 
put; thus at most limited increase of BH mass M, takes place. 

The empirical fact (see Wandel 2002; Vestergaard 2004) 
that a similar M, - a correlation appears to hold also for the 
currently active and faraway BHs is consistent with our adop- 
tion of comparable values for M, energizing both modes of 
nuclear activity: the impulsive quasar feedback effective for 
plasma ejection from groups and galaxies, and the smoother, 
long-lived AGN outputs sufficient to preheat the gas falling 
into poor clusters. We are also consistent with the rough 
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equality of the mass densities in BHs derived from the power- 
ful emissions of the quasars in the optical band (see Marconi 
et al. 2004), and from the weaker and later AGN activity de- 
tected mainly in X-rays (Hasinger 2004; Fabian 2004a). 

On the other hand, strong and rare (i.e., increasingly spo- 
radic) quasar impacts can also explain the scatter of the X-ray 
data widening toward smaller systems as poor groups or mas- 
sive galaxies (Mushotzky 2004). This we trace back to the 
increasing variance in the occurrence of strong quasar events 
or even in their coupling level /, that concur to dynamically 
modulate the plasma ejection Am/m oc fM, and nonlinearly 
affect Lx oc (1 - Am/m) 2 when Am/m approaches 1. As the 
hierarchical clustering proceeds toward clusters, instead, the 
evolution of the quasars cuts down most internal effects; this 
is because the impulsive contributions to AE within a struc- 
ture's dynamical time hardly can keep pace with the increase 
of the equilibrium energy E oc m T v in such late and massive 
systems with deep potential wells. 

We stress that our upper and lower bounds for Lx illustrated 
in Fig. 3 by the SN strip and the quasar line comprise nearly 
all data points, except for a few groups with peculiar features 
currently under scrutiny (see Mushotzky 2004; Osmond & 
Ponman 2004). So we submit that several pieces of data fit 
together when considering both the external preheating from 
AGNs and the internal impact from quasars, with the same 
average values of fM,. We remark that several authors (e.g., 
Ruszkowski & Begelman 2002; Fabian 2004b; Zanni et al. 
2004) have argued the relevance of AGN feedback in explain- 
ing the puzzle posed by the "cool cores" (Molendi & Pizzo- 



lato 2001) at the very centers of many clusters. On the other 
hand, Cavaliere et al. (2002) and Granato et al. (2004) have 
stressed that powerful quasar impacts in massive spheroids 
easily quench star formation and produce precociously red gi- 
ant elliptical galaxies. 

To conclude, we stress that energy feedback from AGNs 
and quasars with an overall coupling around 5 x 10~ 2 to the 
ambient baryons yields agreement with independent observa- 
tions in different frequency bands and over different distance 
scales. Specifically, this paper is focused on the extended X- 
ray emissions and plasma entropy of poor clusters and groups; 
but we have also considered at galactic and subgalactic scales 
the mainly optical correlation of nuclear BH masses vs. host 
velocity dispersions. At the intermediate scales of early mas- 
sive galaxies and in the ^wave/submm band, we have pro- 
posed in Lapi et al. (2003) how to catch quasar impacts in 
the act from resolved Sunyaev-Zel'dovich signals enhanced 
by overpressure in running blast waves. Such transient events 
sweeping plasma outward to lower densities (see Figs. 5 and 
CI) hardly increase the extended X-ray emissions; they in- 
stead specifically correlate with pointlike X-rays from a fully 
active quasar, and/or with strong IR emissions signaling a 
nascent quasar enshrouded by dust. 
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APPENDIX 
A. HYDROSTATIC EQUILIBRIUM 

The hot plasma constituting the ICP/IGP pervades the potential wells of clusters and groups, being in overall virial equilibrium 
with the DM. As the sound crossing time is comparable to, or somewhat shorter than the structure dynamical time, hydrostatic 
equilibrium applies; when the thermal pressure p = nkT//j,is dominant (see Ricker & Sarazin 2001; Inogamov & Sunyaev 2003) 
this yields 

JL^=_^, (Al) 
m p n dr dr 

in terms of plasma number density n, and of the temperature T. The solution of this differential equation requires one boundary 
condition, for example the value n(R) = n% at the virial radius r = R; it also requires an equation of state, i.e., a specific relation 
between n(r) and T(r). 

As to the DM potential <t>(r), we adopt the widely used NFW form (Navarro, Frenk & White 1997) 

m^ g (c) l ^l^, (A2) 
r/r s 

involving the one-dimensional velocity dispersion a of the DM, and the scale r s = R/c. The concentration parameter c w 
5 (M/10 15 Mq)"° 13 slowly increases (Bullock et al. 2001), and the factor g(c) = [ln(l +c)— c/(l +c)]~ l weakly rises from clusters 
to groups. To wit, the smaller, earlier DM halos are more concentrated in terms of the normalized potential <p = $/cr 2 , consistent 
with the tenets of hierarchical structure formation (see Padmanabhan 2003). 
It is useful to recast the hydrostatic equilibrium in terms of the all-important adiabat K = kT/n 2 ^ to obtain 

5 d\nn + dh\K kT v dip _ (A3) 
3 dlnr dlnr Kn 2 ^ dlnr 

recall that kT v = [im p a 2 is the virial temperature. The above is a first-order differential equation of the Euler type for n(r); in 
terms of K(r), it can be formally integrated by standard methods to obtain Eq. (2). 

In the cluster outskirts, Eq. (A3) directly relates the slopes of n(r) and K(r), so that n(r) oc r ~ 3(a+2) / 5 corresponds to K(r) oc r a . 
This is easily derived near rmR, where for the NFW potential dcf>/dlnr w 3 holds; meanwhile the coefficient T v /Kn 2 ^ 3 T v /T2 
is easily recognized to be /3, with values close to | for clusters. 

B. ACCRETION SHOCKS AND LAYERS 

Across any sharp transition like a shock, the conservation laws of mass, total stress, and energy for a plasma with 3 kinetic 
degrees of freedom are written (Landau & Lifshitz 1959) 

«i vi =n 2 v 2 

p i + m p n i v\ = p 2 + nip n 2 \\ 

,3 ,5 _ a 1 ,5 



-m p niv 1 + -pivi = -m p n 2 v 1 +-p 1 V2 



(Bl) 



As is customary, we have indicated with the subscripts 1 and 2 the pre- and postshock variables, respectively; in addition, by v 
we indicate velocities measured in the shock rest frame. 
The previous system of equations leads after some algebra to the temperature jump T%/Ti = 6 under the general form 

e =T6 M ' + rT6j?- (B2) 

in terms of the Mach number M. = (3 /im p v\/5kT\) x l 2 \ this is Eq. (6) of the main text. It is seen that shock heating (8 > 1) 
requires the flow to be supersonic in the shock rest frame, i.e., M. > 1 as expected. The corresponding density jump reads 

n 2 AM 2 

^"^T3 ; (B3) 
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or in terms of 6 (cf. Cavaliere et al. 1999) 

no / 1 \ / / 1 \ 2 1 

(B4) 

V u J u 

This may be used to express the postshock adiabat Ki = kTz/n^ , which leads to Eq. (5) of § 4 in the context of accretion flows. 
The same general equation applies also to the leading shock of an outgoing blast wave, discussed in § 5. 

Focusing now on accretion flows, we consider what happens if the transition occurs across a layer of finite thickness S located 
at r w R. In such a case, the conservation equations above include additional terms due to volume forces or nonplanar geometry; 

e.g., in the momentum equation the gravitational term J^_ s drGM(< r)n(r)/r should be considered. However, these corrections 
are 0(5/ R) relative to the term p\+m p n\v\ when S/R <C 1 applies; in particular, we have checked this to hold for any reasonable 
DM and gas distributions M(< r) and n(r) inserted in the above integral. 

In fact, for the numerous small merging partners with mass ratio M' /M < 5% the diffuse baryonic component is stripped away 
promptly (i.e., within a layer S/R < 10%) from its DM counterpart, and raised to the final temperature and entropy levels (see 
Tormen et al. 2004). These smaller merging partners with cooler T/ undergo prompter and also stronger transitions; so they 
not only are better described by Eqs. (B2), (B3) and (B4), but also dominate the averaging procedure over the distribution of 
progenitor masses M', that are used in the main text to compute the effective value of KziT). For all these reasons, we refer to 
"shocks", but imply that similar results also hold for thin layers. 

It is useful to relate the inflow velocity v\ in the shock rest frame to the infall velocity v\ in the cluster rest frame, to obtain 



V1 = 3 V1 




4 tim p v\ 



(B5) 



here we have assumed the kinetic energy to be small downstream. In turn, the infall velocity v\ is set by the DM gravitational 
potential $2 at the virial radius through energy conservation, which yields 

v, = V"2x$2 ; (B6) 

the fudge parameter x= l-R/Rf expresses the ignorance on the exact position of the radius R/>R at which free fall begins (see 
Voit et al. 2003). The upper bound x ~ 0.7 obtains if Rj is computed by equating the free-fall speed to the Hubble flow. In fact, 
basing on many numerical simulations since Bertschinger (1985) we adopt the value \ ~ 0.37, close to that obtained if the gas 
inflow begins at the turnaround radius during the gravitational collapse of a standard "top-hat" perturbation (see Padmanabhan 
2003). 

We now derive a number of convenient approximations to the temperature and density jumps valid for very strong, moderately 
strong, and weak shocks, that are used in § 4 of the main text. For very strong shocks with kT\ / fim p vj< 1 that occur for 
supersonic accretion onto rich clusters, Eq. (B5) reads v\ ~ 4vi/3; thus Eq. (B2) and Eq. (B4) respectively approximate to 

O^Y-^ ^4; (B7) 

the former expresses the limit where full thermalization of the inflow kinetic energy takes place. We can now use Eq. (B6) with 
$2 ~ -5.1 a 2 (the NFW potential corresponding to a concentration c = 5) to find vi w 2.1 (kT v / ^,m p ) l l 2 ; putting all together yields 
f3 = kT v /9kT\ ~ 3kT v /fim p v 2 w 0.7, a value consistent with the observations of rich clusters. 

For moderately strong shocks with kT\ / yw p v\ < 1 that occur for preheated accretion onto poor clusters, Eq. (B5) gives v\ ~ 
Av\/3 + 5kT\/Aijm p Vi. Now Eqs. (B2) and (B4) approximate to 

a ^mA 3 ^^AisiV (B8) 



3 m 2 m \ 16 9 

relatedly, the parameter [3=T V /9T\ is lowered from the cluster value around 0.7 to 

/J*0.7 (B9) 

For weak shocks with kT\/ \im p v\ » 1 that occur for preheated accretion in small groups, one recovers from Eqs. (B2) and 
(B4) the adiabatic limit 

C. A NEW FAMILY OF SELF-SIMILAR BLAST WAVES 

Hydrodynamic flows are amenable to a self-similar description when their dynamics can be characterized in terms of space 
and time variables, and of a small set of parameters with independent dimensions. Although self-similar solutions are only 
particular descriptions of a hydrodynamic flow, they often accurately yield its actual asymptotic behavior and offer a useful guide 
for understanding its generic features, as discussed by Zel'dovich & Raizer (1967). 
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Self-similar solutions have been systematically investigated by Sedov (1959) for blast waves originated by time-dependent en- 
ergy discharges into a gas with power-law density. Such solutions have then been successfully applied to a variety of astrophysical 
problems, such as the propagation of SN remnants (see Chevalier 1976), and of shocks driven by solar flares (see Parker 1963). 
Most of these treatments consider a gas with negligible initial pressure so that the resulting blast wave is strongly supersonic; 
moreover, they do not include gravitational effects, in particular those due to a dominant DM component. 

In the main text we are interested in strong and weak blasts driven by comparable energy discharges into plasmas in virial 
equilibrium with the DM at temperatures differing by factors 10 2 from galaxies to clusters. To cover the full range, we derive 
here (see also Lapi 2004) a new family of self-similar solutions that not only include a power-law initial density gradient and 
time-dependent energy injection, but also incorporate DM gravity and a finite initial pressure. In Fig. 5 we preliminarily outline 
the unperturbed equilibrium and the perturbed flow. 

C. 1 . The ambient medium 
We consider initial configurations constituted by DM profiles 

f*r) = p R (j) (CI) 

with power-law index 2 < lo < 5/2; the case lo = 2 corresponds to the standard isothermal sphere. The plasma is assumed to 
be in equilibrium within the DM potential well with initial density n(r) oc p(r). Such simple power-law runs will be useful for 
obtaining self-similar blasts, but also constitute fair piecewise approximations of the plasma distributions expressed by Eq. (2). 
The plasma temperature and entropy profiles are given by 

kT _ 2nGp R R" 2 _ u 



[wip (lo — 1)(3— lo) 

= 2nGpm 5 p /3 (p R R^ ^ 
(Lu-l)(3-u)(m/M) 2 / 3r 

(C2) 

relatedly, the polytropic index T = 2(1 - 1/lo) ranges from 1 to 1.2 when lo increases from 2 to 2.5. The total energy of the gas in 
equilibrium is 



£tot« r) = l - 



3 3 ~ W -(4 W -7) 



m(<r)c l s ocr 5 " 2w ; (C3) 



5 5-2lo 

here m(< r) = 4-n m p J r dxx 2 n(x) = 4n (m/M) p R R^ r 3_a, /(3 - to) is the gas mass within r, and the sound speed is c s = 
(5kT /3 pm p ) x l 2 oc r 1 ""/ 2 . The bounds 2 < lo < 2.5 guarantee the total energy of the system to be finite and negative; here- 
after and in the main text we indicate its modulus with E(< r) = -E tot (< r). 

C.2. The blast 

A blast wave sweeps through the plasma as a result of the energy injections AE(t) by a central source. The ensuing unsteady 
gas flow is described by the system of partial differential equations 

d,n + d r (nv)+ — = 
r 

~ ^ ~ 1 * , GM« r) 

o t v + vo r v+ o r p + r = 

m p n r L 

(d t + vd r )-^ = 0, 

(C4) 

supplemented at the leading shock by the Rankine-Hugoniot boundary conditions. The latter may be derived from the general 
expressions Eqs. (Bl) specialized to the case of internal shocks, i.e., taking the preshock flow velocity vi in the shock rest frame 
equal to —v s , opposite to the shock velocity. This yields 

_3 M 2 -l P2 _ 5M 2 -1 » 2 _ 4M 2 
V2 ~4 Vs M 2 pi" 4 n l ~M 2 + 3 1 ( } 

in terms of the Mach number A4 = v s /c s . Note that the condition for shock formation and propagation v s > c s can also be written 
as P2+m p n 2 v\ > p\, i.e., the total stress pushing the shock outward has to exceed the upstream pressure. 

The flow described by the above equations will be self-similar if it can be expressed in terms of the variables r and t and of 
only two more parameters with independent dimensions. One must be the gravitational constant G if gravity effects are to be 
included; as to the other, we take the quantity p R R^ entering the initial state given by Eq. (CI). 

Self-similarity then implies for the law of energy injection 

AE(t) oc G (5 -" )/w (prR^) 5 ^ t 2(5 - 2u>)/u> . (C6) 
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FIG. C 1 . — Distributions of density (top panel), cumulative mass [middle panel) and entropy [bottom panel) within the blast, computed from Eqs. (C9) and 
(C10) for ui = 2, and normalized to their postshock values. Solid lines are for a strong shock with AE/E = 3, dashed lines for an intermediate shock with 
AE/E = 1, and dotted lines for a weak shock with AE/E = 0.3. 



The resulting time-dependencies of the source output turn out to be interesting; the values of the index uj correspond to luminosi- 
ties going from a constant (uj = 2) to a spike (uj = 5/2), the upper range being useful to describe the quasar fading out because of 
its own feedback on the accreting gas. Since E(< R s ) oc R 5 ~ 2uJ oc f 2 * 5 " 2 ")/" cx A£ holds, it is easily seen that AE/E is constant 
during the blast motion, and thus constitutes the key parameter for labeling the solutions. This comes about because the values 
of AE/E sets the Mach number, i.e., the strength of the shock, as we specify below. 

Under self-similarity, Eqs. (C4) are solved along the following lines. First, we use the dimensional parameters of the problem 
to construct the self-similarity (adimensional) variable £ = r/R s (t), where 



R s (t) = 



5-kGprR" u) 2 M 2 



6(w-l)(3-w) 



l/u 



f 2/u 



<xR(Mt) 



2/u 



(C7) 



is the shock radius. It is now seen that the Mach number Ai = v s (t)/c s [R s (t)] is independent of time and position, as is AE/E. 
Note that in our simple model with u = 2 the blast moves out with constant speed, while for to > 2 it decelerates. 
Then we introduce the adimensional quantities V(£)» ^ (£)> an d ^"(0 through 

3/inip r 2 



v(r,t) = -V(0 n(r,t)-. 



1 



T>(0 T(r,t)-. 



-T(Q; 



t ' ' ' m P C,t ~" 5k r 

these enable us to convert the partial differential Eqs. (C4) into a set of ordinary differential equations 



(C8) 



v'+lv- 2 -) % 



= 2-3V 



VI 2 - 

UJ 



V 



3^ (V V 

-5 T (r + v 



6 -T + V 2 - 



5 uj 2 M 2 ^ 



2\ (T_ivr m 
uj T 3 V 



*4 



with boundary conditions at £ = 1 (r = R s ) given by 



V s 



3 M 2 -l 

2 toM 2 



X = 



(5M 2 -l)(M 2 + 3) 
4uj 2 M 4 



24(m/M) (3-uj)(uj-1) 
5^ uj 2 (M 2 + 3) 



(C9) 



(C10) 



As a last step, we have numerically solved the differential system Eqs. (C9) by using a standard Runge-Kutta integrator with 
adjustable time step. For various shock strengths we show in Fig. CI the distributions of density, mass, and entropy in the blast 
for the uj = 2 model. 
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FIG. C2. — The Mach number M (top panel) and the piston position A = R p /R s (bottom panel) as a function of AE/E. The solid line is for ui = 2, the dashed 
line is for ui = 2.25, and the dotted line is for u) = 2.4. 



While the radial mass and entropy distributions differ considerably in the blast and in the initial configuration, the entropy 
distribution n(m) as a function of the integrated plasma mass m(< r) remains unchanged, and reads 

K(m) otm(< r ) (6 -" )/3 (3 " w) . (C 1 1 ) 

The postshock normalization K2 is raised because of dissipation within the shock, after the general Eq. (5) in the main text that 
has been derived in Appendix B. 
The overall energy balance for the blast flow is written 



:R s ) = 4n f ', 
Jo 



AE(l ) — E\ < R,) = 4;v / drr 2 <j ^m p n v 2 + ^p- GM( < r) m[>n K ( C12) 



This equation relates M and AE / E, and the result is shown in the top panel of Fig. C2. Note that for strong shocks AE /E oc M 2 
holds; in particular, such a dependence implies that in the limit M 2 ^> 1 (but in fact already for M 2 > 3) our family of self-similar 
blast waves converge to the standard solutions with negligible gravity and zero initial gas pressure. E.g., for uj = 2 and constant 
source luminosity one has R s oc A4t oc (AE /E) l l 2 t; since E oc R s holds, one recovers the shock motion R s oc L 1//3 f oc A£ ,1//3 f 2/ ' 3 , 
provided by the standard blast wave theory in the strong shock limit. As another example, consider 10 = 2.5 and spiky energy 
liberation after Eq. (C6), for which one has R s oc _A/f 4 / 5 f 4//5 oc (AE /E) 2 I 5 t A / 5 \ since now E = const holds, one recovers the 
standard dependence R s oc AE 2 l 5 t 4 / 5 for strong shock. 

From Fig. CI it is easily seen that the flow is confined within a shell that terminates at the leading shock at R s and begins at a 
trailing "piston", the contact discontinuity located at R p = XR S < R s where the action of the source is transferred to the plasma. 
Self-similarity implies the thickness 1 - A of such a shell to depend only (and inversely) on the shock strength; A is plotted vs. 
AE/E in the bottom panel of Fig. C2. 

Analytic approximations may be derived for the limiting behavior of the adimensional variables V, T>, T in the vicinity of the 
piston. For a given w such limiting behaviors turn out to be independent of the shock strength, and read 

(w-6)/3C7-2u>) 

V 



u) 5 



I- 



1 



T p oc — . (C13) 



Thus at the inner piston the density diverges weakly but the mass vanishes (so the overall effects of radiative cooling are negligi- 
ble), while the temperature goes to zero making up a finite pressure. 

C.3. The shell approximation 

Since the perturbed flow is confined within a shell of constant thickness A between the inner piston and the leading shock, it is 
convenient to represent our solutions by using the shell approximation (see Cavaliere & Messina 1976; Ostriker & McKee 1988). 
Here we improve on the classic treatment by extracting the value of the shell thickness A directly from the exact solution (see 
Fig. C2), in order to obtain results reliable to better than 15%. 
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The equation of motion for the shell is written 



±[m«R s )v 2 ]=4nR 2 s [(p) - Pl ] - GM{ < * s) m{< R s ) ; (C14) 



dt s ' ZJ s 1X1 1 ' IJ 5-2lu R 2 

here (p) is the volume-averaged pressure, which as function of M reads 

M'"»^ llt = . ( C15) 
pi 8 i—uj 5-2w 

For weak shock with M — > 1 we obtain (p)/p\ = 3/(5 -2a;), while in the strong shock limit .M >• 1 our result (p)/p\ — » 
5 7W 2 (8 - 3 w)/8 (3 - w) matches that known for standard, strong blast waves without gravity (e.g., Cavaliere & Messina 1976). 
Integrating Eq. (C14) leads to 

1 , 3 , , 3-uj GM(<R S ) 
AE(t)-E« R s ) =-m« R s ) v z 2 +-(p)V- — m« R s ) , (C16) 

in terms of the shell volume V = 4-irR 3 s (l - A 3 )/3; this is the simplified, shell version of Eq. (C12). 

A relevant quantity is the ratio between the kinetic and the thermal energy of the blast; inserting Eq. (CI 5) into Eq. (CI 6) one 
finds 

^ = I ^fzl ( ci7) 

AE Th 2(8-3w)(l-A 3 ) M 2 

This is easily seen to vanish in the weak shock limit ( Ai — > 1), and to take on values 3/(16-6 w) (1 - A 3 ) for strong shocks with 
M. ^> 1 , as in standard blasts with no gravity. 

Finally, we remark that Eqs. (C7), (C15), (C16), and (C17) constitute a set of handy relations useful in the main text and 
suitable to implement in semianalytic models (SAMs) or hydro+iV-body codes of galaxy formation. 



